Local Statistics of Realizable Vertex Models 

Zhongyang Li* 



Abstract 

We study planar "vertex" models, which are probability measures on edge 
subsets of a planar graph, satisfying certain constraints at each vertex, exam- 
ples including dimer model, and 1-2 model, which we will define. We express 
the local statistics of a large class of vertex models on a finite hexagonal lattice 
as a linear combination of the local statistics of dimers on the corresponding 
Fisher graph, with the help of a generalized holographic algorithm. Using an 
n X n torus to approximate the periodic infinite graph, we give an explicit 
integral formula for the free energy and local statistics for configurations of 
the vertex model on an infinite bi-periodic graph. As an example, we simulate 
the 1-2 model by the technique of Glauber dynamics. 

1 Introduction 

A vertex model is a graph G = {V, E) where we associate to each vertex v & V a 
signature A local configuration at a vertex w is a subset of incident edges of 
V. A configuration of the graph G is an edge subset of G. The signature at a 
vertex f is a function which assigns a nonnegative real number (weight) to each local 
configuration at v. The partition function of the vertex model is the weighted sum 
of configurations X G {0, 1}^, where the weight of a configuration is the product 
of weights of local configurations, obtained by restricting the configuration at each 
vertex. Dimers, loop models, and random tiling models are some special examples 
of vertex models. 

Direct computations of the partition function of a general vertex model usually 
require exponential time. On the other hand, using the Fisher-Kasteleyn-Temperley 
method P1,f5], we can efficiently count the number of perfect matchings (dimer config- 
urations) of a finite planar graph. The idea of generalized holographic reduction 
is to reduce a vertex model on a planar graph to a dimer model on another planar 
graph, essentially by a linear base change, see section 3. For an original version of 
the holographic reduction (Valiant's Algorithm), see [I6] . 

However, not all satisfying assignment problems can be reduced to a perfect 
matching problem ("realized") using the holographic algorithm. We study the re- 
alizability problem of the generalized algorithm for vertex models on the hexagonal 
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lattice and prove that the signature of reahzable models form a submanifold of pos- 
itive codimension of the manifold of all signatures, see Theorem 3.6. An example of 
realizable models is the 1-2 model, which is a signature on the honeycomb lattice, 
only one or two edges allowed to be present in each local configuration, see Figure 
13, 14. The realizability problem of Valiant's algorithm is studied by Cai |2], and 
the realizability problem of uniform 1-2 model (not-all-equal relation), a special 
1-2 model which assigns all the configurations weight 1, under Valiant's algorithm is 
studied by Schwartz and Bruck [Tl]. 

Realizable vertex models may be reduced to dimer models in more than one 
way, that is, using different bases. However, all the dimer models corresponding to 
the same vertex model are shown to be gauge equivalent, i.e. obtained from one 
another by a trivial reweighting. 

One of the simplest vertex configuration models is a graph with the same signature 
at all vertices. Using the singular value decomposition, we prove that such models on 
a hexagonal lattice are realizable if and only if they are realizable under orthogonal 
base change. Moreover, the orthogonal realizability condition takes a very nice form; 
see section 3.2. 

We compute the local statistics of realizable vertex models on a hexagonal lattice 
with the help of the generalized holographic reduction, i.e. for the natural probability 
measure, we compute the probabilities of given configurations at finitely many fixed 
vertices, which are proved to be computable by sums of finitely many Pfaffians, see 
Theorem 5.1 and Theorem 5.2. 

The weak limit of probability measures of the vertex model on finite graphs are of 
considerable interest. Using an nxn torus to approximate the infinite periodic graph, 
we give an explicit integral formula for the probability of a specific local configuration 
at a fixed vertex, see section 6. These results follow from a study of the zeros of the 
characteristic polynomial, or the spectral curve, on the unit torus T^. For a 
more general result about the intersection of the spectral curve with T^, see [13]. For 
example, using our method, we compute the probability that a {001} dimer occurs 
and for uniform 1-2 model, and the probability that a {011} configuration occurs at 
a vertex for critical 1-2 model, see Examples 7.2 and 7.3. 

The main result of this paper can be stated in the following theorems 

Theorem 1.1. For a periodic, realizable, positive-weight vertex model on a hexagonal 
lattice G with period 1x1, assume the corresponding Fisher graph has positive edge 
weights, then the free energy of G is 

F:= hm llog5(G„) = -i^ // \ogP{z,w)- — 

where Gn is the quotient graph Gj (nL^n£) , P{z, w) is the characteristic polynomial. 

Theorem 1.2. Assume the periodic vertex model on hexagonal lattice is realizable 
to the dimer model on a Fisher graph with positive, periodic edge weights, and as- 
sume the entries of the corresponding base change matrices are nonzero. Let A„ be 
the probability measure defined for configurations on toroidal hexagonal lattice Gn. 
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Moreover, for a configuration c at a vertices v 

lim X^icv) = J][n«;,J|Pf(X-^)v(,^.) 

where 

-f^o^^((«> yi), {v, X2, y2)) = 



1=1 



47r2 



^xi-x2yjyi-y2 



f2 



P{z, w) iz iw 



dj are local dimer configurations on the gadget of the Fisher graph corresponding to 
V. V{dj) is the set of vertices involved in the configuration dj. 
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2 Background 

2.1 Vertex Models 

Let {0, 1}*^ denote the set of all binary sequences of length k. A vertex model is a 
graph G — (V, E) where we associate to each vertex v eV a, function 

: {0, lY'^^""^ R+ 

r^ is called the signature of the vertex model at vertex v. We give a linear ordering 
on the edges adjacent to v, and we fix such an ordering around each vertex once 
and for all. This way the binary sequences of length deg(t') are in one-to-one cor- 
respondence with the local configurations at v. Each edge corresponds to a digit; if 
that edge is included in the configuration, the corresponding digit is 1, otherwise the 
corresponding digit is 0. Hence we can also consider r^, as a column vector indexed 
by local configurations at v: 





( r„(0. 


.00) \ 




r.(0. 


.01) 


rv = 


r.(0. 


.10) 




^ r,(l. 


■.11) ) 



Example 2.1 (signature of the vertex model at a vertex). Assume we have a degree-2 
vertex with signature 



( r.m \ 

r.(Ol) 
r.(lO) 

V ^^(11) / 



7 

\^ J 



that means we give weights to the four different local configurations as in Figure 1: 
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Figure 1: Relation at a Vertex 

Assume G is a finite graph. We define a probability measure with sample space 
the set of all configurations, Vt = {0, l}'^'. The probability of a specific configuration 
7^ is 

m = ^llr.,{n). (1) 

The product is over all vertices. r„(7^) is the weight of the local configuration ob- 
tained by restricting TZ to the vertex v, and S" is a normalizing constant called the 
partition function for vertex models, defined to be 

The sum is over all possible configurations of G. 

Now we consider a vertex model on a Z^-periodic planar graph G. By this we 
mean that G is embedded in the plane so that translations in act by signature- 
preserving isomorphisms of G. Examples of such graphs are the square and Fisher 
lattices, as shown in Figure 7. Let G„ be the quotient of G by the action of nZ^. It is 
a finite graph embedded into a torus. Let Sn be the partition function of the vertex 
model on G„. The free energy of the infinite periodic vertex model G is defined to 
be 

F := lim ^logS'„ 

n->oo n 

2.2 Perfect Matching 

For more information, see [6]. A perfect matching, or a dimer cover, of a graph 
is a collection of edges with the property that each vertex is incident to exactly one 
edge. A graph is bipartite if the vertices can be 2-colored, that is, colored black 
and white so that black vertices are adjacent only to white vertices and vice versa. 

To a weighted finite graph G = (V, E, W), the weight W : E ^ R+ is a function 
from the set of edges to positive real numbers. We define a probability measure, called 
the Boltzmann measure fi with sample space the set of dimer covers. Namely, for 
a dimer cover D 

^{D) = I n ^(^) 

e&D 

where the product is over all edges present in Z), and Z is a normalizing constant 
called the partition function for dimer models, defined to be 

Z = J2llW{e), 

D eSD 
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the sum over all dimer configurations of G. 

If we change the weight function W by multiplying the edge weights of all edges 
incident to a single vertex v by the same constant, the probability measure defined 
above does not change. So we define two weight functions W to be gauge equiv- 
alent if one can be obtained from the other by a sequence of such multiplications. 

The key objects used to obtain explicit expressions for the dimer model are 
Kasteleyn matrices. They are weighted, oriented adjacency matrices of the graph 
G defined as follows. A clockwise-odd orientation of G is an orientation of the 
edges such that for each face (except the infinite face) an odd number of edges point- 
ing along it when traversed clockwise. For a planar graph, such an orientation always 
exists [5]. The Kasteleyn matrix corresponding to such a graph is a |V^(G')| x |V^(G')| 
skew-symmetric matrix K defined by 



It is known [H |5l [151 E] that for a planar graph with a clock-wise odd orientation, 
the partition function of dimers satisfies 



Now let G be a Z^-periodic planar graph. Let Gn be a quotient graph of G, as 
defined before. Let '~^x,n{ly,n) be a path in the dual graph of Gn winding once around 
the torus horizontally(vertically). Let Eh{Ev) be the set of edges crossed by jxily)- 
We give a crossing orientation for the toroidal graph G„ as follows. We orient all 
the edges of G„ except for those in Eh U Ey- This is possible since no other edges 
are crossing. Then we orient the edges of E^ as if Ey did not exist. Again this 
is possible since G — Ey is planar. To complete the orientation, we also orient the 
edges of Ey as if Eh did not exist. 

For 6,T E {0, 1}, let -R'^'^ be the Kasteleyn matrix Kn in which the weights of 
edges in Eh are multiplied by (—1)^, and those in Ey are multiplied by (—1)"^. It is 
proved in [15] that the partition function Zn of the graph Gn is 



Let Em = {ei = UiVi, ...,em = UmVm} be a subset of edges of G„. Kenyon [7] 
proved that the probability of these edges occurring in a dimer configuration of Gn 
with respect to the Boltzmann measure P„ is 



where E^ = V{Gn) \ {ui, vi, Um, Vm}, and {K^'^)e^ is the submatrix of fC^"" whose 
lines and columns are indexed by E!^. 

The asymptotic behavior of Zn when n is large is an interesting subject. One 
important concept is the partition function per fundamental domain, which is defined 



u,v 




Z = Vdet K. 
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to be 

lim \ log Zn. 

n— s>oo n 

Let Ki be a Kasteleyn matrix for the graph Gi. Given any parameters z,w, 
we construct a matrix K{z,w) as follows. Let •jx,!, ly,i be the paths introduced 
as above. Multiply K^^v by z if Pfaffian orientation on that edge is from u io v, 
otherwise multiply K^^^ by ^, and similarly for w on 7^^. Define the characteristic 
polynomial P{z,w) = det K{z,w). The spectral curve is defined to be the locus 
P{z,w) = 0. 

Gauge equivalent dimer weights give the same spectral curve. That is because 
after Gauge transformation, the determinant multiplies by a nonzero constant, and 
the locus of P{z,w) does not change. 

A formula for enlarging the fundamental domain is proved in [3l [8]. Let Pn{z, w) 
be the characteristic polynomial of with period 1x1, and Pi{z,w) be the char- 
acteristic polynomial of Gi, then 

Pn{z,w)= II II Pi{u,v) 

u"-=z v"=w 

2.3 Matchgates, Matchgrids 

A matchgate F is a planar local graph including a set X of external vertices, i.e. 
vertices located along the boundary of the local graph. The external vertices are 
ordered anti-clockwise on the boundary. F is called an odd(even) matchgate if it has 
an odd(even) number of vertices. 

The signature of the matchgate is a vector indexed by subsets of external 
vertices, {0, Ij'^L For a subset Xq C X, the entry of the signature at Xq is the 
partition function of dimer configurations on a subgraph of the matchgate. The 
subgraph is obtained from the matchgate by removing all the external vertices in 
Xo. 

Example 2.2 (signature of a matchgate). Assume we have a matchgate with external 
vertices 1, 2, 3, and edge weights as illustrated in the following figure: 
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Figure 2: Matchgate 
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then the signature of the matchgate is 



/ 000 \ 




/ 


ax + by + cz 


\ 


001 











010 











oil 






X 




100 











101 






y 




110 






z 








V 





/ 



A matchgrid M is a weighted planar graph consisting of a collection of match- 
gates and connecting edges. Each connecting edge has weight 1 and joins an external 
vertex of a matchgate with an external vertex of another matchgate, so that every 
external vertex is incident to exactly one connecting edge. 



3 Generalized Holographic Reduction 

In this section we introduce a generalized holographic algorithm to compute the par- 
tition function of the vertex model on a finite planar graph in terms of the partition 
function for perfect matchings on a matchgrid. The idea is using a matchgate to 
replace each vertex, and perform a change of basis, such that after the base change 
process, the signature of a vertex becomes the signature of the corresponding match- 
gate. We describe the algorithm in detail as follows. 

For a finite graph G, we associate to each oriented edge e a 2-dimensional vector 
space Vg. To the edge with the reversed orientation, the associated vector space is 
the dual space, i.e. V^e — K*- G:ive a set basis {/^ , } for each Vg, satisfying 

fUfi) = s,,. 

Let w be a vertex with incident edges e;^, e/^, oriented away from v. The signature 
of the vertex model at a vertex v, r^, can be considered as an element in — 
Vei^ ® ... ® Vei^. Hence r„ has representations under bases F — {fe:fe}eeE and 
S = 6g}eg£; as follows 

= ry{ci, ■ ■ ■ cijbeil ® • • ■ ® bl\l (2) 

where B are the set of standard bases for each 




Ci. G {0, 1}. Q-^ ■ ■ ■ Qj. are binary sequences of length k. From the definition of the 
signature of vertex models, obviously r^(Q^ • • • Q^) is the weight of the configuration 
• • • cii^ at vertex v. 



7 



We construct a matchgrid M as follows. We replace each vertex f by a matchgate 
Pt,, such that the number of external vertices of is the same as the degree of f, 
and the edges of the vertex model graph G become connecting edges joining different 
matchgates in the matchgrid M. Examples of such replacements are illustrated in 
the following Figure. 




Figure 3: degree-3 vertex and matchgate 




Figure 4: degree-4 vertex and matchgate 
If the signature m„ of the matchgate satisfies 

= ^ r^,/(Qi • • ■ ci^)lfei\ ® ■ ■ ■ ® h%l , (4) 

that is, the representation of under bases B is the same as the representation of 
r„ under bases F, then we have the following theorem: 

Theorem 3.1. Under the above base change process, the partition function of the 
vertex model of G is equal to the partition function of the dimer model of M. 

Proof. There is a natural mapping $ from ^y^yW^ to C induced by ^eeE4>e, where 
(pe is the natural pairing from Ve^V* to C. Note that in ^v^vW^, each Ve and V^e 
appear exactly once. Since the representation of niy under bases B is the same as 
the representation of r„ under bases F, we have 

$(®^,ev"^^,) = ^{®vevrv)- (5) 

([5]) follows from the fact that each 0e : K ® K* — ^ C is independent of bases as long 
as we choose the dual basis for the dual vector space. However, the left side of ([5]) is 
exactly the partition function of the dimer model of M, while the right side of (j5]) is 
exactly the partition function of the vertex model of G. □ 
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Define the base change matrix at edge e, Tg = ( /g /g ). The base change 
matrix at vertex f , T^, acting on Wy by multiplication, is defined to be 

®{g|g is incident to v, and oriented away from i)}-^e- 

In order for a vertex model problem to be reduced to dimer model problem, one 
sufficient condition is that at each vertex, the signature of the vertex model Vy under 
the bases F is the same as the signature of the matchgate m„ under the standard 
bases. Namely, 

TyTriy = Ty. (6) 

([6]) follows directly from (|3]) and (j4]). 

Example 3.2. Consider the graph in Figure 3 with standard dimer signature 



a 




Figure 5: 



rw = rb = ( rooo ^"001 '^oio '^oii '^loo '^loi '^iio '^iii f 

= (01101000)* 
Define the base change matrix on edges a, b, c to be 



Ta = T, = Tg 



1 

1 



By definition, Ty^ = Ta®Ti,®Tc. Note that Ty = iTl) = Tb After the base change, 
we have the standard loop signature 

'^W Ty; • Tyj Xf, " T }j 

= (00010110)* 

For instance, after the base change, the dimer configuration 001 with only c-edge 
occupied becomes 

Ty, ■ (bo ® 60 ® bi) 

= (00000010), 

which is the configuration 110, the loop configuration with a-edge and b-edge occupied. 



9 



However, since the number of vertices in a matchgate is either even or odd, 
at a vertex of degree d, the signature of a matchgate must be a 2*^"^ dimensional 
subspace of , with those 2'^~^ entries being 0. These are the entries correspond 
to the partition function of dimer configurations on a subgraph of the matchgate 
with an odd number of vertices, see example 2.1. This is the parity constraint. 
As a result, by dimension count we can see that it is not possible to use holographic 
algorithm to reduce all vertex models into dimer models. To characterize the special 
class of vertex models applicable to holographic reduction, we introduce the following 
definition. 

Definition 3.3. A network of relations on a finite graph is realizable, if there exists 
a system of bases reducing the model to the set of perfect matchings of a matchgrid. 

Definition 3.4. A network of relations is bipartite realizable, if it is realizable 
and the corresponding matchgrid is a bipartite graph. 

Remark. The above generalizes Valiant's algorithm [16] in the sense that our basis 
can be different from edge to edge. As a consequence, our approach results in an 
enlargement of the dimension of realizable submanifold, which will be shown in the 
next section. 

3.1 Realizability 

We are interested in periodic vertex models on the honeycomb lattice with period 
n X n. The quotient graph can be embedded on a torus T'^ = x S^. We classify 
all the edges into a-type, 6-type and c-type according to their direction, and assume 
6-type and c-type edges have the same direction as the two basic homology cycles 
(1, 0) and (0, 1) of torus, respectively, see Figure 6. 



Assume the vertex model is realizable, then each corresponding matchgate may 
have either an even or an odd number of vertices. By enlarging the fundamental 
domain, we can always assume that there are an even number of matchgates with 
an even number of vertices. Then by permuting rows of matrices on a finite number 
of edges, we can always have all the matchgates having an odd number of vertices. 
For example, assume we have a pair of adjacent matchgates with matrix on the 




(0,1) 



■'(1,0) 



Figure 6: Periodic Honeycomb Lattice 
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connecting edge e, Tg = ( J° ) (Tg here is actually the inverse of the base change 

ti 
to 



ti 

matrix defined before ), as illustrated in Figure 7. Then if we assume 





Figure 7: Permutation of Basis Vectors 1 

obtained by permutating two basis vectors of Tg, we actually interchange the roles of 
and 1 at the corresponding digit of the binary sequences as indices of signatures of 
the matchgate. Namely, assume m„ = T„r„, and m„ = T^r^, where Ty = Ta®Tb®Tc. 
Then for any binary sequence C1C2C3, the entry my^cic2C3} is the same as ^d{(i-ci)c2C3}) 
according to equation ([6]). If originally we have an even matchgate at f, by parity 
constraint, the 001,010,100,111 entries of are zero. After the permutation of 
basis vectors, fh^ will have 101,110,000,011 entries to be zero. Hence has to 
be an odd matchgate, as illustrate in Figure 8. By finitely many times of such 





Figure 8: Permutation of Basis Vectors 2 

permutations, we can move two even matchgates to adjacent positions. Then we 
permutate the basis on the connecting edge, we can decrease the number of even 
matchgates by 2. Since we have an even number of even matchgates in total, if we 
repeat the process, in the end, all matchgates will be odd. 

Assumption 3.5. From now on, we make the following assumption: 

• All entries of base change matrices are nonzero. 

• All entries of the matchgate signature are nonzero. 

Since the honeycomb lattice is a bipartite graph, we can color all the vertices in 
black and white such that black vertices are adjacent only to white vertices and vice 
versa. Assume vertex signatures at black and white vertices are as follows: 

ly^J —— { rf^^3 /y»^J /-v»^J ryi^J | 7 I 

' w — \ "''OOO -''OOl '''OlO '''Oil "''lOO "''101 "''110 "''111 j \' J 

ij / ij ij ij ij ij ij ij \^ {Q\ 

~ { Vooo Vooi Voio Voii Vioo Vioi Vuo Viii ) \°) 



11 



where < i < n,l < j < n is the row and column index of white and 

black vertices. Assume 000 = 1, then entries of signatures are indexed from 1 to 
8. Associate to an a-edge, 6-edge and c-edge incident to a white or black vertices, a 
basis 

({i,j,k) {i,j,k) \ 

nff'"^ j - 1, ^1 Pi ) 

I (i,j,k) (i,j,k) \ 

rp(hhk) _ / P02 \ — I „{id,k) (i,j,k) 

~ I (i,j,k) (i,j,k) \ — \ ^2 P2 

^12 Pl2 J ^ 

(i,j,k) {i,j,k) \ 

7^(ij,fc) _ I '^OS P03 \ — I y^ihiM) ^{id,k) 

^13 Pit / ^ 



where k = 1 for white vertices and = for black vertices. 

By the realizability equation (Q, and the fact that all matchgates are odd, the 
1st, 4th, 6th, 7th entries of the matchgate signatures m„ are 0, so we have a system 
of 4 algebraic equations at each vertex v. For example, at each black vertex, we have 

m,,,,o = T^''''®n'^''®Tl:'^''-rl' (9) 
and the fact that the 1st entry of rriij^ is gives the following equation 

noino2no3yi + ^01^02^031/2 + ^0lP02^03l/3 + ^0lP02P03l/4 

+Poino2no3y5 + Poino2Po3y6 + PoiPo2no3y7 + P01P02P03I/8 = 

That the 4th, 6th, 7th entries of friij^ are give three other similar equations. 
Similarly at each white vertex, we have 

m,,,i = [(T^'^'^ ® Tl'^^' ® T^'^'^Y]-' . (10) 

the same process gives a system of 4 equations at the white vertex. 

Let of = "(f - 1) = -ffjoT; ^ ^ {0, l},m G {1, 2, 3}. Then the equations we get 

m 

are linear with respect to a^i, af\, 6q^, which can be solved explicitly. 

(11) 
(12) 

(13) 

(14) 
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f^J . y^J 
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■ v^^ 
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. fi 






K^i ■ 




lij 
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. fi 
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where 




Since aom 7^ ciim and bom 7^ ^im, after clearing denominators and some reducing, 
for any equations (|TT1) -( |T4|) are equivalent to 

27/27/8 - 27/67/4 + (613 + &03)(l/lZ/8 + ?/2|/7 " 1/51/4 " VeVs) + 2(7/l7/7 - 2/5Z/3)&13&03 = 0(16) 

-2x1X7 + 2x5X3 + (ai3 + ao3)(a;2a;7 + ^iXg - X4X5 - X6X3) + 2{xqX4 - X2Xs)aisao3 = 0(17) 

21/81/3 - 27/41/7 + (&12 + hQ2){,ym + 1/31/6 - 1/41/5 " I/7I/2) + 26i2&02 (l/el/l " I/2I/5) = 0(18) 
2X5X2 - 2xiX6 + (xiXs - X2X7 - X5X4 + X3X6)(ai2 + 002) + 2(X7X4 - X3X8)ao2a-12 = 0(19) 

If we solve ao2, 012, &02, &12 explicitly, a similar process yields 
2x3X2 — 2x1X4 + (X5X4 + xiXs — X2X7 — X3X6)(aoi + ail) + 2(x6X7 — X5X8)aoiciii = 0(20) 

2i/6i/7 - 2i/5i/8 + (1/3I/6 + 1/21/7 - ym - i/4i/5)(&oi + bii) + 2(1/21/3 - yiy4)boibii = 0(21) 

We get two equations per edge, one involving the a-variables, the other involving 
the b- variables. But a-variables and b-variables are actually the same thing for each 
single edge. From equation (fT6ll . (ITTl) . we can solve basis entries all ~ ^« 3~^! from 
equation (fT8i) . (fT9ll . we can solve basis entries a]'l = b]'^^'-' . Finally from (|TTl) -( fT^ . 

we can solve a;_i and the only constraint left will be a]f = b^'l, which is two 
polynomial equations with respect to relation signature at each a-edge. Together 
with Theorem 2 in appendix, we have the following theorem 

Theorem 3.6. Under assumption 3.5, the realizable signatures on the n x n peri- 
odic honeycomb lattice form a lAn"^ dimensional submanifold of the IGn"^ dimensional 
manifold of all positive signatures; the bipartite realizable signatures on the nxn peri- 
odic honeycomb lattice form a Vlr? dimensional submanifold of the 16n^ dimensional 
manifold of all positive signatures. 

For the exact realizability condition, see the appendix. 

Under assumption 3.5, the weight {111} is non-vanishing at each matchgate. 
Since the probability measure will not change if all entries of the signature at one 
vertex are multiplied by a constant, we can assume at each matchgate, the weight of 
{111} is 1. Therefore, we have 

Theorem 3.7. A realizable vertex model on a finite hexagonal lattice G can always be 
transformed to dimers on M , a Fisher graph as shown in Figure 9, with the partition 
function of the vertex model of G equal to the partition function of the dimer model 
of M , up to multiplication of a constant. 
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Figure 9: Matchgrid with 3x3 Fundamental Domain 

It is possible to construct a matchgrid with different weights to produce the same 
vertex model. Since the holographic reduction is an invertible process, by which 
we mean that because the base change matrices are nonsingular, we can always 
achieve the matchgate signature from the vertex signature and vice versa, there is 
an equivalence relation on dimer models producing the same vertex model. 

Definition 3.8. A vertex model is holographically equivalent to a dimer model on 
a matchgrid, if it can be reduced to the dimer model on the matchgrid using the 
holographic algorithm, in such a way that partition function of the vertex model cor- 
responds to perfect matchings of the matchgrid. Two matchgrids are holographically 
equivalent, if they are holographically equivalent to the same vertex model. 

Proposition 3.9. Holographically equivalent matchgrids give rise to gauge equivalent 
dimer models. Therefore, they have the same probability measure. 

Proof. See the appendix. □ 
3.2 Orthogonal Realizability 

Definition 3.10. A vertex model is orthogonally realizable if it is realizable by an 
orthonormal base change matrix on each edge. 

Consider a single vertex. To the incident edges of the vertex, associate an ma- 
trices Ui,U2jUs- Without loss of generality, we can assume Ui,U2,U^ G 5*0(2). In 
fact, if detUi = —1, we multiply the first row of Ui by —1. The new signature of 
the matchgate will be multiplied by —1. This change does not violate the parity 
constraint. 



14 



Assume Ui = i sma^ \ Assume U = ®^_i?7j, then U G 5*0(8), each 

y — sm Oi cos Oi J 

term of which is product of three of cosa^jSinaj. Moreover, the eigenvalues are 
g(±ai±02±a3)v^^ Using trigonometric identities, each entry is a hnear combination 
of cos(±ai ± a2 ± as) and sin(±ai ± 0:2 ± as). If we further define 



9 = { 9i 92 93 94 95 96 97 98 )' 
h = [0 h2 h-s /i5 hs)' 
7 = tti + 0:2 — as; = "1 + "3 — "2; = a2 + as — ai 
ji = 9i + 96 + 97 - 91, h = 91 +96 + 97 - 9i\ 
is = 9i+ 9a + 97 - 96] ji = 91+94 + 96 - 97] 
J5 = 93 + 95 + 98- 92] J6 = 92 + 95 + 98- 93] 
37 = 92 + 93 + 98 - 95] J8 = 92 + 93 + 95 - 98 



P = j2 cos(v3) + Jt sin(v3); Q = js cos{4j) + jg sin(?/'); R = J4 cos(7) + sin(7) 

K = ji cos{(p + + 7) - js sin(v9 + + 7) 

H = j-j cos(v9) - j2 sin((y9); L = cos(V') - js sin(V'); M = jg cos(7) - sin(7) 

iV = js cos(v3 + ip + 'y) + ji sm{(p + ■?/; + 7) 

then if Ug = h, we have 

= P = Q = i? = ir (22) 

Ah2 = H + L-M + N (23) 

Ahs = H-L + M + N (24) 

Ah = N + M + L-H (25) 

Ah8 = H + L + M-N (26) 

By (I22D, we have tanv? = -4^,tan^ = -4a,tan7 = -^,tan(v9 + ^ + 7) = 4i. If 
we define 

u + V 



t{u, v) 
tt{a, b, c, d) 



1 — uv 

a + b + c + d — abc — abd — acd — bed 



1 ~ ab — ac — ad — be — bd — cd + abed 
Then the following theorem holds: 

Theorem 3.11. A vertex model on a periodie honey eomh lattiee with period n x n 
is orthogonally realizable if and only if its signatures satisfy the following system of 
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equations 



tt( 



z. 



ijk 
7 



ijk 



Z. 



ijk 



ijk 

Z( . 



ijk ' ijk ' ijk ' ijk 



+ ( ^6 



6 



^2 



£4 \ _ 
TJtr; - 

5 

f2 

7 

j4 



Zr 



t 



^8 

^6 _ 



z-i 1 
24 • 



2 



f 



^2 



6-, 



where z^ 



ijl 



ijl 
■ Xi, Z2 



f( f7 ^4 

^2 

a:3 + a;5 + a;8 



X4 + X6 + X7 

ijl 

Xq — Xj, Zq 



ijl 
X2, Z^ 



ijl 
■ X3, Z^ 



X2 + X5 + Xs 

Xi + X4 + Xj — Xq, Zj^'^ = Xi + X4 + 
X2+X3 + X5—XS, and the same relation for z^-^^ and yi, ...,ys- x's and 



y's are defined by ^ and 

It is trivial to verify these equations in any given situation. 

Definition 3.12. A vertex model on a periodic honeycomb lattice with period n x n 
is positively orthogonally realizable if it is orthogonal realizable and for each vertex 
{i,j,k), there exists angles ^p^^^ ^ip"^^^ ^--^^^^ , such that 



smcp 



sin-?/' 



sin7 



sin(— 7 — {p — ip) 







zf 








ijk 






{zf? 


ijk 

z( 





+ iz.. 



ijk\2 



Proposition 3.13. // a vertex model on a bi-periodic hexagonal lattice is orthogo- 
nally realizable, then the corresponding dimer configuration has positive edge weights. 

Proof. Under the assumption that the vertex model have nonnegative signature, we 
have 

^1 ^ ^4 + ^6 + Zj] ^4 ^ ^1 + ^6 + Zy] Zq <^ Zi -\- Z4 -\- Zj] Zi <^ Z4 -\- Zq -\- Z-j 
Z2<Z^ + Z^ + Zs] Z3<Z2 + Zr, + Zq] Z^ < Z2 + Z3 + Zq] Zg < Z2 + Z3 + Z^ 

at all vertices. Since at least three of ^1,^4,^0, 2:7 are nonnegative, similarly for 
2^2, -23, 2^5, 2^8, if we take absolute value for all Zi, the above inequalities also hold. 
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therefore 



zi + zi + \ zi + zi + \ zi + 2;8" - ./^2 + > 



^4 + ^5 + V^l + ^8 + V^2 + ^7 - \/^3 + ^6 > 



zi + 2^8^ + ^2^6^ + zi + ^^2^ + ^7 - + ^5 > 

Under the assumption that both the vertex model are positively orthogonal realiz- 
able, the left side of the above inequalities are exactly edge weights of corresponding 
matchgates. □ 

Theorem 3.14. // all matchgates have the same signature, for a generic choice 
of signature, a vertex model on a hexagonal lattice is realizable if and only if it is 
orthogonal realizable. 

Proof. Obviously orthogonal realizability implies realizability, we only need to show 
that realizability implies orthogonal realizability. 

Assume r — [ Xi X2 x^ x^ x^ Xq Xj Xg ) is a realizable signature for the 
vertex model. By lemma 4.1, we can assume the corresponding bases on all edges 
have real entries. Consider the singular value decomposition for the base change 
matrix on each edge. Without loss of generality, assume 

where Ui,Vi G 0(2), and Aj is a nonnegative real number. Assume 
V = (<H)i=iVi)r* ^ {vi V2 V3 V4 V5 ve v-j Vg Y 

Then we can consider ((H)f=i ^ g \.^)'^ (<^f=i ^ q ?^ ^ to be signatures of 

black and white vertices that are orthogonal realizable. Then by the conditions of 
orthogonal realizability, we have 



Vl + X1X3VQ + X1X2VJ - A2A3V4 _ ^1 + aTa;^6 + 

X3V2 + X2V3 + A1A2A3V8 - X1V5 ■^V2 + ^^^3 + xJ-.x.^ Vs - j^V^ 

111 

vi + AiA3t;6 - X1X2V7 + A2A3t;4 _ ^1 + aTa^^s - 171^^7 + x^x;^^. 



-X3V2 + X2V3 + AiA2A3T;8 + A1W5 -^^^2 + ^^^3 + aT^^S + Xr^5 

Vl - X1X3VQ + AiA2f7 + A2A3t;4 _ ^1 ~ l^^e + ArA;^7 + A^^4 

A3t;2 - A2t;3 + AiA2A3f8 + Ai^s -^V2 - jr^V3 + J^J^^Vq + 

-Vl + AiAgfg + X1X2V7 + A2A3t;4 -Vl + xkl'^a + A^^7 + xhr^^ 



X3V2 + A2V3 - A1A2A3V8 + X1V5 ^V2 + J^V3 - x^J^Vs + J-^ 



V5 



(27) 
(28) 
(29) 
(30) 
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If we transform fractions into polynomials, and consider (|27l) + (|28|)-(|29|)-(l30|). (|27l) + (l29|)-(l28|)-(l30|), 
and (l27!) + (l30])-((28])-(l29]), then 

(Ai + l)(Ai - 1){V4VS + V1V5 + V3V7 + V2V(i) = 

(A2 + 1)(A2 - l){v3Vi + v^ve + V1V2 + VjVs) = 

(A3 + 1)(A3 - 1){V4V2 + V-jV^ + V'iVi + fgfe) = 

Then for a generic choice of signature, we must have Ai = A2 = A3 = 1. □ 



4 Characteristic Polynomial 

Assume all entries of the vertex signatures are strictly positive. In this section we 
prove some interesting properties of the characteristic polynomial. 

Lemma 4.1. For realizable vertex model with positive signature and period 1x1, 
there exists a realization of base change over GL2(M), (i. e. one can take base 
change matrices to be real), with the property that, at each edge, < 0. 

Proof. For 1x1 fundamental domain, aik = bik, V /, k. By (lTT]) -( JT4l) . we have 



p 


■t 


s 


■t 


r 


■ u 


Q 


■ u 




■ w 


s 


■ w 


r 


■ V 


Q 


■ V 


(1- 


t 


r ■ 


' t 


s ■ 


u 


_ p- 


u 




w 


r ■ 


w 


s ■ 


V 


p- 


V 



From which we derive 



{jp ■ t){r ■ v) + {p ■ w){r ■ u) — {r ■ t){p ■ v) — {r ■ w){p ■ u) = (31) 

(s ■ t){q ■ v) + {s ■ w){q ■ u) — {q ■ t){s ■ v) — {q ■ w){s ■ u) = (32) 

{p ■ t){q ■ v) + {p ■ w){q ■ u) — {q ■ t){j> ■ v) — {q ■ w){j> ■ u) = (33) 

(s ■ t)(r ■ v) + {s ■ w){r ■ u) — {r ■ t){s ■ t>) — (r ■ w){s ■ u) = 0. (34) 



Since each base change matrix is invertible, we have aoi 7^ an. Plugging (JTSll into 
(!3Tl) -( |34l) . and factor out (ao2 — o-u), or (aos — ctis), we obtain that 093, 013 are two 
roots of quadratic polynomial 

{xeVb + Vy^s + yiX2 + y'iXi)z^ + (-2/7X7 - 2/3X3 + y2X2 + y^x^^ + x^ya - x^y^ + ygXg - yiXi)z 

-yex^ - X32/4 - y2Xi - xjys = 0, 

and ao2, 012 are two roots of quadratic polynomial 

{x^yi + yaxs + y2X4 + y->xy)z^ + {y^x^ - x^yQ - x^y^ + 2/7X7 - y2X2 + ysx^ + 2/4X4 - yiXi)z 

-ysX6 - yjx-, - X12/3 - y4X2 = 

Under the assumption that Xi,yj are positive, these polynomials have real roots, since 
the products of two roots are always negative. Then the lemma follows. □ 
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Figure 10: Weighted 1x1 Fundamental Domain 

By definition, the characteristic polynomial P{z, w) is the determinant of an 
6x6 matrix K{z, w) whose rows and columns are indexed by the 6 vertices in a 1 x 1 
fundamental domain, see Figure 8: 



P{z, w) = det 








Cl 


-bi 


-1 








\ 




-Cl 





ai 





_ 1 

z 









bi 


— ai 











w 






1 






















z 





-C2 





0.2 




v 








w 


b2 


-a2 





/ 



X X 2 uj 

= (z + -){ab - c) + (w + — )(ac - 6) + (— + — )(6c - a) + + 6^ + 1 + 
z w z 

where 

a = 0x02 b = 6162 c = C1C2 (35) 

and we consider the non-degenerate case, which means ac — b ^ 0, be — a 0, 
ab — c ^ 0. We have the following lemma: 

Lemma 4.2. For realizable vertex model with period 1x1, let a, b, c, d be the product 
of dimer weights [SB^) after holographic reduction, under assumption 3.4, we have the 
following inequalities: 

{a + b — c — d){a + c — b — d){a + d — b — c) >0 

a + b + c + d>0. 

Proof. For generic choice of vertex signature, we can assume quotients of basis entries 
Q-ij = — ~ are finite. Apply the realizability equation (|6]) to a 1 x 1 quotient graph, we 
have 8 equations for each black vertex, and 8 equations for each white vertex. At each 
black vertex, 4 equations has on the right side, and 4 equations with ai, 61, ci, di on 
the right side. We divide the 8 equations into 4 groups, each of which has 1 equation 
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with on the right, 1 equation with a non- vanishing edge weight on the right. We 
take the difference of the two equations in each group, and we get 4 new equations 
with a non- vanishing edge weight on the right. We perform the same procedure for 
each white vertex. Then we multiply and add those equations correspondingly, we 
obtain 

a + d-b- (-ao3 + ai3)(y32^4 + XsUv + Xey^ + yiX2) 
a + c-b- (ao2 - ai2){yiX3 + y2X4 + y&xs + y->X7) 
a + b- c- d^ (aoi - aii)(yiX5 + y^xs + y2X& + y^x^) 

Moreover, since holographic reduction leaves the partition function invariant, 

a + 6 + c + d = xxyx + X2y2 + a^^aZ/s + 2:4^4 + x^y^ + x^y^ + xyyr + xsy% 

where the left side is the partition function of dimers of the 1x1 quotient graph, and 
the right side is the partition function of the vertex model, see Figure 10. According 
to lemma 5.1, aofe and aijt are two real numbers with opposite sign. Without loss 
of generality, we can assume and Aqs f-rc both positive. At the black vertices, 
the {000} entry of matchgate signature corresponds to partition functions of perfect 
matchings of the generator with all output vertices kept. Since there are odd number 
of vertices, no perfect matching exists, therefore 

aoi = ■ ■ ■ < 

002^031/1 + a02l/2 + ^031/3 + 

we have 

-(ooi - aii)(ao2 - cii2)(ao3 - 013) > 

The lemma follows from the assumption that entries of relation signatures are posi- 
tive. □ 

The expression of the characteristic polynomial shows that P{z, w) is a smooth 
function on = = 1, |w| = 1}. We are interested in the intersection of 

spectral curve P{z,w) — with unit torus T^, because it has implications on the 
convergence rate of correlations. Theorem 4.4 describes the behavior of P{z, w) on 
T^. Before stating the theorem, we mention an elementary lemma: 

Lemma 4.3. Let /(0) = Asm(f) + Bcoscf) + C, where A,B,C are real. If C > 
and A^ + B"^ -C^ < 0, then f{(f)) > for any e R. 

Theorem 4.4. Under assumption 3.4, either the spectral curve after holographic 
reduction is disjoint from T^, or intersects at a single real node, that is, for some 
{zo,Wo) = (±1,±1), 

P{z,w) = a{z - zof + /3{z - Zo){w - Wq) + -f{w - Wof + 

where — Aa^ < 0. 
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Proof. Let 



[2{ac — 6) + 2 (6c — a) cos 9] cos + 2(6c — a) sin 6 sin ( 



+ 2{ab - c) cos 9 + + + + 1. 



Consider Q{9, (j)) as a trigonometric polynomial with respect to 0, 

(5(6', 0) = ^ sin + S cos + C 

where 

A = 2(ac - 6) + 2(6c - a) cos 61 
5 = 2{hc-a) sin^ 

C = + 6^ + + 1 + 2(a6 - c) cos^ > 0. 

Define 

g{cos9) = C'' - - 

= A{ah-cf cos^ ^ + 4(a6 + c) (a^ + 62-0^-1) cos ^ 
- 4(6c - a)2 - 4(ac - hf + (a^ + ^2 + ^2 + 1)2, 

then gi is a quadratic polynomial with respect to cos 9 with discriminant 

A = 64a6c(a + 6 - c - l)(a + 6 + c + l)(a + 1 - 6 - c)(a + c - 6 - 1). 

The minimal value of g{t) is attained at 

(a6 + c)(a2 + 62_c2-i) 



in = 



where 



^(^o) 



2{ah - c)2 
A 



16 (a6 - c)2 
If a6c > 0, 

^(cos 9) = [2{ah + c) cos ^ + + 6^ - - 1]^ + 16a6c sin^ ^ > 

therefore Q{9, 0) > 0, and Q{9, (f)) — only if sin ^ = sin = 0. 

If abc < 0, lemma 4.2 implies that A < 0, then 5i(cos 9) > 0, therefore Q{9, 0) > 0. 
If abc — 0, we have I ^1^1 ~ 1. Moreover, 

' I ab—c I ' 

(a^+6^-c2-l)^-4(a6-c)2 = {a+l+b+c){a-l-b+c){a+l-b-c){a-l+b-c) > 0, 
which implies 



ab + c 



a' + b'-c'-l 



2{ab - c) 



> 1 



ab — c 

then g{cos9) > 0, and Q{9,(f)) > 0. 

Therefore the only possible intersection of the spectral curve with is real point. 
It is trivial to check that ^L,, = 0, f^L,, , ,x =0, and Hessian\P(z,w)] is 

positive definite wherever a real intersection of P{z, w) = and exists, then the 
theorem follows. □ 
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5 Local Statistics 



5.1 Configuration at Single Vertex 

We are interested in the probability of a specified configuration at a single vertex. 
Consider a realizable vertex model on a planar finite hexagonal lattice, whose parti- 
tion function is S. Assume both the vertex with specified signature and its neighbors 
lie in the interior of the graph. We work out an explicit example here; other cases 
are very similar. If only the configuration {011} is allowed at a black vertex v, we 
get a new vertex model, assume partition function is S'ooi ■ Then the probability that 
the local configuration {001} appear at v is 

Pr({001},t;) = ^. 

Hence the local statistics problem reduces to the problem of how to compute S'ooi 
efficiently. We will show that S'ooi can also be computed by the dimer technique with 
the help of the holographic algorithm. 

Let us denote the adjacent vertices of v by Va,Vb,Vc, according to the direction of 
the connecting edges, as illustrated in Figure 11. 



Va<^ 




Figure 11: Configuration Oil at a Black Vertex 

Let ry, Ta-i Tb, Tc denote the signatures of relations at vertices Va, vj,, Vc, then 
we have 
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Here ~k means the entry at the position can be arbitrary. This is because we only 
allow the configuration {011} at v, therefore the only configurations which actually 
affect the partition function of the vertex model will be those who do not occupy 
the a-edge of Va, and occupy both the b-edge for and the c-edge for v^. We will 
split each of ra,rb,rc into 2 parts, namely, ri = r/(0) + r;(l), I = a,b,c, below. By 
definition, the partition function of the vertex model with signature r,ra,rb,rc is 
a sum of 8 terms, each of which is the partition function of a vertex model with 
signature ry,ra{i),rb{j),rc{k) i,j,k e {0, 1}. That is, 

ij,fee{o,i} 

For each S{ry,ra{i),rt{j),rc{k)}, we give new bases on edges adjacent to v, such that 
^{rv,ra{i),rb{j),rc{k)} bccomc the partition function of certain local dimer configurations. 
Namely 

„ base change „ /q'7\ 

'~'{rv,ra{i),rb{j},rc{k}} — ^{mvii,j,k),ma(i},nn,(j),mc(k)} ["^ ' ) 

where Z{my(^ij^k),ma{i),mb{j),mcik)} is the partition function of the dimer model on a 
matchgrid with signature m^^i, j, k) , rUaii) , mb{j) , mc{k) . Assume the original match- 
gates Ua, Ub, Uc have weights 

mi — [ a2 &2 4 0^2 )*, l — a,b,c 

see Figure 12. We require 




Figure 12: Weighted Matchgate 
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Notice that for any / G {a,b,c}, m/(0) corresponds to configurations with an 
unoccupied /-edge; and mi{l) corresponds to configurations with an occupied Z-edge. 
We are trying to determine the new base change matrices on incident edges of v, for 
each rv,ra{i),rb{j),rc{k), such that fl36l) . fl37|l . fl38l) are satisfied simultaneously. 

Assume on the adjacent edges of Va, ff,, fc we have original base change matrix 

On an a-type edge adjacent to v, we consider 2 possible base change matrices 
qo _ f ^01 \ qi _ f 

V %i (111 J V %i 

On a b-type edge adjacent to v, we consider 2 possible base change matrices 



On a c-type edge adjacent to f , we consider 2 possible base change matrices 



03 "ns 

' " V %°3 ; ^3-^0 ql. 

For signatures ra{i),rb{j),rc{k), we choose basis S{, S^, on edges a, b, c. By 
the realizability condition at vertex v, we have 

m^{i,j, k) = {SI (g) 5^ (g) S^y ■ = e{ijk}Wijk, 

where 

Wijk = m\^q]2<ikzyii 

and ejjjfc} is the 8 dimensional vector with entry 1 at the position labeled by the 
binary sequence {ijk}, and elsewhere. Obviously the choice of new base change 
matrices, namely, the position of zeros in the new base change matrices, results in the 
single configuration {ijk} at the matchgate corresponding to v. Now the problem is 
to determine the entries of the base change matrices satisfying the equations. 
According to realizability conditions at vertices Va,Vb,Vc, 

Ta = T.lM^) = ^LiSl®T^^T^-m,{{) (38) 



n = T.]^M3) = ^UTi®Si®Tlmb{3) (39) 

= Tl^^r,{k) = Tl^{ri®T^®Sl-m,{k). (40) 

For the original basis, we have the realizability condition as follows: 

ri = Ti®T^®T^-mi. (41) 
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Substituting r^, r^, Tc by (|4T]) in ( l40l) we have the following system of linear equations 
with respect to entries of 5* 
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Under the assumption that original base change matrices are invertible, we have 
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Therefore, we only need to choose the nonzero entries of the new base change matrices 
to be the same as the original ones. 

For the other seven configurations at vertex v, we can use the same method to 
achieve a similar result. By splitting each one of r^, r^, into 2 parts, we express the 
partition function of the satisfying assignments as the sum of 8 terms, each one of 
which is the partition function of the relations with signature r^, ra{i), rb{j),rc{k), 
k G {0, 1}. We apply the holographic reduction to each part separately, and we 
derive that the partition function of the relations with signature ra{i), rb{j), rdk), 
i,j,k G {0,1} is equal to the partition function of the dimer configurations on 
the corresponding matchgrids with signature m^j^i, j, k),ma{i),mb{j),mc{k), which 
corresponds to the local configuration {ijk}. This may not be a dimer config- 
uration; to make it a dimer configuration, let us divide those 8 terms into two 
groups: one consists of all {i,j,k} satisfying i + j + k = 0( mod 2), namely 
{000}, {Oil}, {101}, {110}; the other consists of all {i, j, k} satisfying i + j + k = 1{ 
mod 2). The sum of partition functions in the first group is equal to the partition 
function of dimer configurations in a matchgrid with weights illustrated in the left 

-; the sum of partition functions corre- 



graph of Figure 11, where t = - , 

o -f o ' uiQii +u;ioi +10110 ■ 

sponding to {001}, {010}, {100}, {ill} is equal to the partition function of dimer 
configurations in a matchgrid with weights as illustrated in the right graph of Figure 
11, up to a multiplication constant wm- 

However, in the left graph we change the parity of the total number of vertices, 
as a result, there is no dimer. Therefore we have the following theorem 

Theorem 5.1. Assume we have a vertex model Hq on a finite hexagonal lattice, 
which is holographic equivalent to the dimer model on a Fisher graph J^q. If at a 
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Figure 13: Even and Odd Matchgates 



fixed interior vertex v ofHo, only one local configuration is allowed, then the partition 
function of the new vertex model Hi can be computed by the partition function of the 

dimer model on a Fisher graph T\, multiplied by a non-zero constant. JFq and T\ 
have the same edge weights except for the matchgate corresponding to v. 

5.2 Configuration at Finitely Many Vertices 

We are interested in the probability that one specified configuration is allowed at 
each of the given finitely many vertices for a realizable vertex model on a hexagonal 
lattice. For simphcity, assume all of them are black. Let Vq be the set of all black 
vertices with specified configuration. To compute the probability, we first give a cri- 
terion to construct the new base change matrices on incident edges of Vq according 
to whether the edge is present in the dimer configuration and in the vertex model 

I is the original base change matrix on the edge, 
Po Pi J 



Base Change Matrices 


Presence in Vertex Configuration 


Presence in Dimer 


\ Po Pi J 


No 


Yes 


/ no \ 
\ Po Pi J 


No 


No 


f no ni \ 
\ pi J 


Yes 


Yes 


[no Hi \ 
\P0 J 


Yes 


No 



whether the edges incident to vertices in Vq are present in the configuration of the 
vertex model is known, given all the specified configurations at Vq. As before, we 
are going to split the signature of each adjacent white vertex into several parts, and 
we apply the holographic reduction to each parts separately. Our expectation is that 
after the reduction process, each part will be equivalent to the partition function of a 



26 



single local dimer configuration. The new base change matrix on each edge is chosen 
according to whether we want the edge to be present in the dimer configuration or 
not after the reduction. As before, all the non-vanishing entries of the new bases 
change matrices are equal to the original ones; we will see how it works below. 

Consider an arbitrary white vertex w G r(V"Q''), the neighbors of Vq. Assume r„,, is 
the vertex signature at w, and is the original dimer signature at the corresponding 
matchgate. Assume 



Let D{w) denote the number of adjacent vertices of w with specified configuration. 
We classify T{Vq) according to D. 

If D{w) = 1, we split = r^(l) + r^(2). Without loss of generality, assume the 
left-digit edge(a-type edge, horizontal edge) connects w with b G Vq, and the edge 
wb is present in the configuration of the vertex model. Assume 



( noi 


nn \ 




' no2 


ni2 \ 




V Poi 


Pn J 




K P02 


Pl2 J 





no3 ni3 

P03 Pl3 



are original base change matrices on the edges adjacent to w. According to the table, 
the presence of wb in the vertex model configuration implies two possible choices of 
the new basis on the horizontal edge, namely, 

^ ( "-01 nil \ = ( 

^ \Poi J 1 V pii 

while on the two other incident edges of w we have the original bases T2, T3. Assume 



m. 



m. 



w 



(1) = ( b^ ) 

(2) = ( Cy, dy,) 



Notice that the nonzero entries m^(l) have indices {001} and {010}, which corre- 
sponds to configurations without the a-edge present; the non-vanishing entries of 
m^(2) have indices {100} and {111}, which corresponds to configurations with the 
a-edge present. Choose 

Then we can check 

r«;(l) + ryj{2) = • = 

If D{w) = 2, we split = '}2t=i''^w{i)- Without loss of generality, assume the 
middle-digit (b- type edge) and right- digit (c-type edge) edges connects w to &2,&3 £ 
Vq, wb2 is present while wbs is not present in the configuration of the vertex model. 
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According to our criteria listed in the table, on wb2 we have two different base change 
matrices 

On W63 we have two different base change matrices 

Si = f ^03 \ , ^ / n,3 \ 

V ^^03 Pl3 / V ^^03 Pl3 / 

while on the left-digit edge, we keep original basis 71, the original basis on the middle 
and right digit edge are T2, T3, as in the D{w) — 1 case. Assume 

m^(l) = (0a^ 000000) 
m^(2) = ( 6,,, ) 
m^(3) = ( c^, ) 
m^(4) = (0000000 d^) 

Obviously, m^(4) corresponds to both edge being present, mw{^) corresponds to 
neither being present, 171^(2) corresponds to b-cdgc being present and c-edge not, 
and niyu^l) corresponds to b-edge not present and c-edge present. Again according 
to the table, choose 

r^(l) = Ti (8) ® Si ■ m^(l) 
r^(2) = Ti (g) S^ ® • m^(2) 
r^(3) = Ti (8) ® • m^(3) 

we can check 

4 

i=l 

If D{w) = 3, we split rw = Ylt=i''^w{'i), and assume myj{i),i = 1,...,4 as in 
D{w) = 2. Without loss of generality, assume Sl, Sf as in D{w) = 1, S2, 3^,31, S^ 
as in D{w) — 2; that is, we have vertex-model signature {110} at w. Choose 

r^(l) = Sl S^ Si ■ m^(l) 
r^(2) = Sl ® ^3^ ■ m^(2) 
r»(3) = 5i2®52^®5^m^(3) 
r^(4) = S^ ® ^2^ ® Si ■ m^(4) 

again we can check = X]i=i D(ty) = 2. 
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For all the other local configurations, the same technique works, and we will have 
a similar result. Assume Vq — {vi, Vp}, r(V^'') = {wi, Wk}, then 



where f^^ is the specified configuration at the vertex Vg. The first equality follows 
from the definition of the partition function of satisfying assignments. When we 
compute the tensor product of relation signatures of all black (white) vertices, we get 
a vector of dimension 2l-^('^)l, where |-E'(G')| is the number of edges in the planar finite 
graph G. This vector can be indexed by binary sequences of length \E(G)\. Each 
binary sequence corresponds to a configuration, and the entry there is the product 
of weights of the local configurations obtained from the configuration restricted to 
each black (white) vertices. Obviously the inner product of the vector at black 
vertices and the vector at white vertices are exactly the partition function of the 
satisfying assignments. The second equality follows from the multi-linearity of the 
tensor product. For 1 < j < fc, if D(wj) — 1, ij G {1,2}; if D{w.j) = 2,3, ij G 
{1,2,3,4}. For each summand, we choose a basis on incident edges of vertices in 
Vq according to whether the edge is present in the dimer configuration and relation 
configuration, and keep the original bases on all the other edges, as described above. 
This is realizable because for each b e Vq, whether its incident edges are occupied 
by dimers are completely determined in each part of the sum. In other words, each 
term on the right side of the second equality corresponds to an configuration on all 
the edges incident to vertices in Vq, which is a local dimer configuration at each odd 
matchgate corresponding to white vertices w's with D{w) = 3, To make them be 
dimer configurations at each black matchgate, we divide those configurations into 
groups according to the following criterion: two configurations arc in the same group 
if and only if the parity of the number of occupied incident edges at each black vertex 
m is the same. If zero or two incident edges are occupied, we construct an even 
matchgate with modified weights at the black vertex; if one or three incident edges 
are occupied, we construct an odd matchgate with modified weights. Notice that 
the modified weights depend only on the local configuration on the incident edges 
of the vertex. Since \Vq\ = p, we have 2^ different constructions in total, but only 
2P~^ of them admit dimer cover, each of which has even number of even matchgates. 
Therefore we have 

Theorem 5.2. Assume we have a realizable vertex model "Ho on a finite hexagonal 
lattice, which is holographic equivalent to the dimer model on a Fisher graph Fq. 
Let Vq be a subset of black vertices. If for each v e Vq, we only allow one local 
configuration, this way we obtain a new vertex model T-Li. The partition function 
of Hi is equal to the sum of partition functions 0/2^" dimer models on matchgrids 
J-*!, ...,J^2P-^, each of which has an even number of even matchgates. J^i{l < i < 2^~^) 
are the same J-'q, except for the matchgates corresponding to vertices in Vq. 





ii,...,ik 
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5.3 Ising Model and Vertex Model 



Consider the Ising model on a finite Kagome lattice, embedded into a torus, as 
illustrated in the following figure. 




Figure 14: Kagome Lattice and Honeycomb Lattice 

The associated honeycomb lattice is illustrated in the Figure by the dashed line. 
Each vertex of the Kagome lattice corresponds to an edge of the honeycomb lattice; 
hence each Ising spin configuration on the Kagome lattice corresponds to an edge 
subset of the honeycomb lattice. If the spin is "+", then the corresponding edge is 
included in the subset; otherwise the edge is not included. Assume the bonds of the 
Kagome lattice have interactions Ji, J2, J3, as illustrated in Figure 14. We define a 
vertex model on the honeycomb lattice with signature at all vertices as follows. 
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This way the probability measure of the Ising model on the Kagome lattice is equiv- 
alent to the probability measure of the vertex model on the honeycomb lattice. It is 
trivial to check that this vertex model is orthogonally realizable. 

6 Asymptotic Behavior 

In this section, we prove the main theorems concerning the asymptotic behavior of 
realizable vertex models on the periodic hexagonal lattice, as stated in the introduc- 
tion. Consider an infinite periodic graph G, with period 1x1, see Figure 9. Our 
technique to deal with such a graph is to consider a graph Gn with 1 x 1 fundamen- 
tal domains, embed G„ into a torus, and consider the limit when n — > 00, see Page 
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5. Our first theorem is about the free energy of the infinite periodic hexagonal lattice. 



Proof of Theorem 1.1 Assume M„ is the corresponding matchgrid with respect to 
Gn, and Pn{z,w) is the characteristic polynomial. Obviously M„ is also a quotient 
graph of a periodic infinite graph modulo a subgraph of generated by (n, 0) and 
(0,n). The corresponding Kasteleyn matrices here are defined given the orientation 
of Figure 9. For even n, the crossing orientation can be obtained from the orientation 
of Figure 9 by reversing all the 2;-edges and w-edges. By Theorem 3.6, M„ is a Fisher 
graph, and the partition of the vertex model can be expressed as follows according 
to the principle of holographic reduction: 

S{Gn) = Z{Mn) = ^1 - Pfif„(l, 1) + FfKnih -1) + Pfi^„(-1, 1) + Pfi^n(-1, -1)| 

Thus 

max |PfP„(M,f)| < Z(M„) < 2 max |PfP„(M,t;)| 

u,ve{-i,i} «,'"e{-i,i} 

On the other hand, according to the formula of enlarging fundamental domains, 
^log max \PiKn{u,v)\= max log\P{z,w)\ 

By Theorem 4.4, either P{z,w) has no zeros on T^, or it has a single real node, 
in which case any sample point in maXu^^g{_i^i} ^^z^=u^^n=jjhg\P{z,w)\ is at 
least — from the real node, for some constant C > 0. Therefore 

lim ^S{Gn) = lim max ^ i \og\P{z,w)\ 

Iff dz dw 

log Piz, w) — - — 



□ 

For each a measure A„ is defined as in ([1]). Assume at a fixed vertices f , we 
only allow configuration c^. Let /i^ be the Boltzmann measure of dimer configura- 
tions on Mn- Let be the matchgrid corresponding to the vertex model which 
only allows at f, as described in Theorem 5.1. Let m^, be the matchgate of Mn-, 
corresponding to f , and dj be a local dimer configuration at m^,, Wd^ be product of 
weights of matchgate edges included in dj, and Vd^ be the set of external vertices of 
matchgates iriy occupied by dimer configuration dj. In our graph, every matchgate 
has 3 external vertices. Our second theorem is about the asymptotic behavior of the 
measure A,,. 
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Proof of Theorem 1.2 According to Theorem 5.1 

Z(M„) _ ^ Z{Mr,{d,)) 

^ Z(M„ \ Vjd,)) 



where the sum is over all local dimer configurations dj. Since the number of local 
configurations is finite, it suffices to consider lim„_>oo ^'^^^ziM^)^'''^ • Note that M„ 
differs from M„ only on edge weights of m^, hence M„ \ V{dj) and M„ \ V{dj) are 
the same. 

Given dj, let Wn^^ = Mn \ V{dj) be the subgraph of M„ by removing all vertices 
occupied by the conffguration {dj), as well as their incident edges. Then 

§^ = ^1 - ^ml\w^4^)) + Pf 

Wi{Kl^-\Wn4,)) + Pf(/C''"'(w^n,.,))| 

First of all, let us assume P{z,w) has no zeros on T^. According to the formula of 
enlarging fundamental domains, for any (6*, r) G { — 1, 1}, Pf{K^'^) ^ 0. Then 

|Pf(X^)^.| = |Pf(i^^)^i||Pf(ii'^)| 

In [21 E], it was proved that given two vertices and {y,X2-,y2) 

Since P{z,w) has no zeros on T^, we have 



ra-s-oo 4:71 J J 12 P[Z,W) IZ IW 

As — 7- oo, each entry of {K^)~^ is convergent, so is the Pfaffian of a finite order 
sub-matrix {K^)y^^^^ y and we have 



Z(M„ \ V{d,)) 



If P(-2, ti?) has a real node on T^, without loss of generality, we can assume the 
real node is (1, 1). It was proved in [Ij that if P{z,w) = has a node at (1, 1), then 
for any fixed finite subset E 



n-s>oo 2Zn 



and 



n— >-cx) 2Zn 



If we take E = V{dij, dpj), our theorem follows. □ 
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7 Simulation of 1-2 Model 



Assume we have 1-2 model with signature = c b a a b c ) at all 
vertices. In other words, at each vertex, either one or two edges are allowed to 
be present in a local configuration. By Section 3.2, the signature is orthogonally 

realizable with base change matrix T = ( '"^^ i_ A. | on all edges. 

\^ - sm ^ cos ^ y 

We define a discrete-time, time- homogeneous Markov chain ^fftt with state space 
the set of all configurations of the 1-2 model. For an n x n honeycomb lattice 
embedded into a torus, the state space is finite, and let us denote it by {ik}^^i- Let 
T{ik) be the set of configurations that can be obtained from configuration ik by adding 
or deleting a single edge. Assume p, q, r are edge variables, namely p,q,r ^ {a, b, c}, 
and {p,q,r} = {a,b,c}. Define Tp^^q{ik) be the set of configurations of 1-2 model 
which can be obtained from i^. by adding a single g-edge uv; before adding uv, only 
a p-edge is present at both u and v. Define Tp_g{ik) be the set of configurations 
of 1-2 model which can be obtained from by deleting a single g-edge uv; after 
deleting uv, only a p-edge is present at both u and v. Define To{ik) = T{ik) \ 
{^p,q&a,b,c,pjtq^p,+q{'ik) ^p,q&a,b,c,pj^q ^p-q{^k)}- Dcfiuc cutrics of transition matrix for 
DJlt as follows: 

if ii e To{ik) 
if ii e Tp^+g{ik) and r>p 
if ii G Tp^+g{ik) and r <p 
if ii e Tp and r <p 

if ii e Tp and r > p 

if ii = ik 
else 

Obviously, OJtt is aperiodic. For more information on Markov chain, see [HI [12] . 
Moreover, we have 

Proposition 7.1. DJlt is irreducible. 

Proof. By definition, we only need to prove that any two configurations communicate 
to each other. We claim that any two dimer configuration can obtained from each 
other by finite steps. In fact, the symmetric difference of any two dimer configurations 
is a union of finitely many loops. Obviously one dimer configuration can be obtained 
from any other dimer configuration by first adding finitely many edges to achieve 
their union, then deleting alternating edges of loops. Notice that we have a 1-2 
configuration at each step. Hence we only need to prove that any configuration of 
1-2 model can reach a dimer by finite steps, each of which is adding or deleting one 
single edge. 

Let us start with an arbitrary configuration of 1-2 model. There are 3 types of 
connected local configurations: loops with even number of edges; zigzag paths with 
odd number of edges; zigzag paths with even number of edges. For the first and 
second types, we can always achieve dimers by deleting alternating edges. Hence 
we can assume that all the zigzag paths with even number of edges are of length 2, 



P(ii\ik) 



1-E 



3n2 
1 

3^ 
1 p2 

3r? 

jj6r{»fc) 





P{ 
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and all the other connected local configurations are dimers. There are two types of 
length- 2 paths. One has vertices black-white-black (BWB), and the other has vertices 
white-black- white (WBW). For each fixed configuration, the number of BWB paths 
is the same as the number of WBW paths; otherwise the complement graph cannot 
be covered by dimers. Consider an arbitrary WBW path in a fixed configuration, as 
illustrated in Figure 12. 



One can easily check that no matter what the configuration is, it communicates 
with a configuration satisfying one of the following two conditions: 1. the number 
of length-2 edges is decreased by 2; 2. it can be moved to a WBW configurations 
containing Bi, for all 1 < i < 6. Hence if the number of length-2 paths does not 
decrease, one can transverse all black vertices without meeting a BWB path, because 
our graph is finite. However, this is impossible because a BWB path always exists 
as long as a WBW path exists. □ 

For an irreducible, aperiodic Markov chain DJlt with transition matrix P and 
stationary distribution tt, let a^o be an arbitrary initial distribution. Then 



Therefore, in order to sample a configuration, we can approximately sample accord- 
ing to the distribution XqP^, with large A^. To that end, first we choose a fixed 
dimer configuration with probability 1 as the initial distribution Xq, then we ran- 
domly change the configuration by adding or deleting a single edge according to 
the conditional probability specified by the transition matrix P. If neither adding 
nor deleting U1U2 ends up with a satisfying configuration, we just keep the previous 
configuration; else we get a new configuration by adding or deleting uiu^- Then we 
repeat the process for N steps. This way we get a sample for distribution XqP^ , 
if is sufficiently large, this is approximately a sample for distribution tt, which is 
exactly the distribution given by 1-2 model. 

Example 7.2. (Uniform 1-2 Model) Consider the 1-2 model with a — h — c — 1. 
After the holographic reduction, the signature of each matchgate is 




Figure 15: WBW Path 



lim xqP"^ = TT. 
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It is gauge equivalent to a positive-weight dimer model on Fisher graph, whose spectral 
curve does not intersect T^. We are interested in the probability that a {001} dimer 
occurs, that is, at a pair of adjacent vertices Vi,V2 connected by an a-edge, only the 
configuration {001} is allowed. By the technique described in section 5 and section 
6, the partition function of the configurations in which a dimer is present at V1V2 is 
equal to the sum of two partition functions, one corresponds to both vi and are 
replaced by an odd matchgate with signatures 



( ^ ^ ); 



the other corresponds to both vi and V2 are replaced by an even matchgate with 
signatures 

v{2 A A _a/2 a _V2 V2 n 
4^'-' 4" 44" 

Both of them have the same partition function as a graph with positive weights. More- 
over, if we give a base change matrix ( ^ ^ ] on V1V2 edge, and apply holographic 



reduction, we see that the two graphs with positive weights are holographic equiva- 
lent. Hence it suffices to consider the one with a pair of odd matchgates with modified 
weights. Then 

Pr{a dimer is present at V1V2) = 

ISZ 

where Zqqi is the partition function of dimer configurations with weight 1 on triangles 
corresponding to vi and V2 and weight | on all the other triangles, and Z is the 
partition function with weight | on all the triangles. Moreover 

2;^^^ — _^ooiooi _|_ ^001111 _|_ ^^010010 ^010100 _|_ ^100100 _|_ ^100010 _|_ ^111001 _|_ ^111111 
_ ^^001001 _|_ 2_^ooiiii _|_ 2_^oiooio _|_ 2_^oioioo _|_ ^111111 

where Z'^^^''^^^ is the partition function of dimer configurations with fixed configuration 
ijk,ijk on triangles corresponding to vi and V2. The second equality follows from 
symmetry. Meanwhile 

12 2 2 

j!^ vOOlOOl _|_ vOOllll _|_ f_ vOlOOlO , ;^ vOlOlOO _|_ viiiiii _ v 

9 3 9 9 

then 

1 / g^'Ooiooi 4_2'0°i^i^ -^Q^Qimio ■j^g^'Oioioo 
Pr(a dimer is present at V-1V2) — — ( 1 H 1 — 1 — h 



18 V 9Z 3Z 9Z 9Z 

1/4 4 Ifi 

(1 + -^\{K-J)2,\ + 9|Pf(i^-^)2356| + y |(i^-^)l3| 
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where 

, 3 ff w{zw + 1 + z — 9w) dz dw fT^-i\ 

i 00)23= JJ^^ 2(;22 + ^2) + 2{z + w) + 2zw{z + w)- 21zwJ^nJ " ^ °° '''' 



13 



w{w + z + z^ — 9zw) dz dw 

IStt^ J Jj2 2{z'^ + w"^) + 2{z + w) + 2zw{z + w) — 21zw iz iw 

y 00)25- JJ^^ 2(^2 + ^2) + 2(z + w) + 2zw(z + w)- 2lzwl^l^ ~^ °° 

1 _ 1 /"/ -z-zw-w + 9 ^^^'^_(T^-i\ 

y 00 J26 - JJ^^ 2(^2 + ^2-) + 2{z + w) + 2zw{z + w)- 2lzwl^^ " ^ ~ '^'^ 

T/ie entries of the inverse matrix can he expressed as elliptic functions. The proba- 
bility that a {001} dimer occurs is approximately 6%. A sample of the uniform 1-2 
model is illustrated in Figure 16. 

Example 7.3. (Critical 1-2 Model) Consider the 1-2 model with a — A,b — c — 1. 
After the holographic reduction, it has the same partition function as a positive-weight 
dimer model on Fisher graph, whose spectral curve has a single real node on T^. The 
probability of the configuration {011} 

Pr({01l}) = ^ = Im^'U + |(i^-^)23| + |(i^-^)l3|) + |Pf(i^-^)l23456|) 



Zqii is the partition function on a Fisher graph with weights 1,1,1 on one triangle, 
221 

3' 3' 3 



and weights 1,1,^ on all the other triangles. Z is the partition function with weights 



|, |, I on all the triangles. By symmetry |(i^^^)i2| = |(i^go^)23|- Then 

Pr{{on}) = ^ + ^|(i^-^)i2| + ^|(i^-%| 



1 1 
+ 



3 67r2 
1 



37r 



w{4:wz + 4 + z — 9w) dz dw 

'j2 —7{w'^ + z'^) + 32wz{w -\- z) -\- 32{w + z) — llAwz iz iw 
z{9wz — Az — w"^ — Aw) dz dw 

fz —7(^2 + z"^) + 2>2wz{w -\- z) -\- 32(w + z) — llAwz iz iw 



1 1 

+ 



3 127r 



18 16^2 - 85^ + 29 

H =dz 



\Z\: 



1 

+ 3^ 



.1 32z - 7 z{32z - 7) ^(32^ - 1)^ Az'' - 11 z ^ A 

9 A 22^2-35^ + 8 

H , dw 



,1 {32w - 7) w{32w - 7) w{32w - 7)VAw^ - 17w + 4 

23 25 4 65 44 

= arctan arctan pa 39%. 

48 1127r 3 3 367r 117 

For \f- we choose a branch with positive real part. This integral can be evaluated 
explicitly because the graph has critical edge weights. A sample of the critical 1-2 
model is illustrated in Figure 17. 

Question: How large should be as a function of the size of the graph? 
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Figure 16: Sample M Uniform 1-2 Model 



Figure 17: Sample of Critical 1-2 Model 

A Realizability Conditions 

In this section we give explicit realizability conditions for periodic honeycomb lattice 
embedded on a n x n domain. Let us consider an arbitrary vertex x with adjacent 
vertices w. Assume relations have following signatures: 

= ( xi ... xs y ry^{ yi ... ys )* 
{ zi ... zs y ry,= { Wi ... ^8 )* 

According to the realizability equation (1), after a lengthy computation, we have 

Theorem A.l. A periodic network of relation onnxn honeycomb lattice is relizable 
if and only if at each vertex x with adjacent vertices y, z, w connected by a, b, c edge, 
respectively, the signatures satisfy the following equation: 

J2 XijkYiZjWk = 
ij,fee{o,i,2} 

where 

Yo = ym - V2yz Yi = y^yg + y^y^ - 2/22/7 - 2/32/6 >2 = 2/52/8 - 2/62/7 

Zq — ZiZq — Z2Z5 Zi — ZiZs + ZzZq — Z4Z5 — ZiZ-j Z2 — Z^Zs — z^z^ 

Wo — WiWj — Wi — WiWs + W2W-J — W5W4 — W3WQ W2 — W2W8 — WqW4 
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and 



Xqoo = xI{x3Xq + ,T2.T7 + X/^X., - .Xi.Xg) - 2X1X2X3X5 

Xool = xl{xQXi - X2Xs) + xl{xiX'r - X3X5) 

X002 = xl{x2Xj - XiXs - X3X6 - 0:4X5) + 2X1X2X4X6 

^010 ~ X^(xiX6 — X2X5) -|- X^(x7X4 — XsXg) 

Xoii = XiX2(x4X7 - X3X8) + X3X4(xiX6 - X2X5) 

X012 = X2(X4X7 - XsXs) + X4(XiX6 - X2X5) 

X020 = X3(X3X6 - X2X7 - X4X5 - XiXs) + 2X1X4X3X7 

X021 = X4(xiX7 - X3X5) + X3(X6X4 - X2X8) 

X022 = X4(X2X7 + XiXs + X3X6 - X4X5) - 2X2X3X4X8 

Xioo = a;5(xiX4 - X2X3) + xKxqXj - XgXg) 

Xioi = XiX2(x6X7 - X5X8) + X5X6(xiX4 - X2X3) 

X102 = xl{XiXi - X2X3) + X2(X6X7 - X5X8) 

Xlio = XiX3(x7X6 - X5X8) + X5X7(xiX4 - X2X3) 

^111 = X1X4X6X7 - X2X3X5X8 

X112 = X4X8(xiX6 - X2X5) + X2X6(X4X7 - X3X8) 

X120 = X7(xiX4 - X2X3) + X3(X6X7 - X5X8) 

X121 = X3X7(X4X6 - X2X8) + X4X8(xiX7 - X3X5) 

X122 = X8(xiX4 - X2X3) + X4(X6X7 - X5X8) 

X200 = X5(X4X5 - XiXs - X3X6 - X2X7) + 2X1X5X6X7 

X201 = X5(X6X4 - X2X8) + X6(xiX7 - X3X5) 

X202 = 2^6(2^1^8 + X2Xy + X4X5 - X3X6) - 2X5X5X2X8 

X210 = xI{x4Xy - X3X8) + X7(xiX6 - X2X5) 

X211 = X6X8(xiX7 - X3X5) + X5X7(X4X6 - X2X8) 

X212 = Xg(xiX6 - X2X5) + xl{x4X7 - X3X8) 

X220 = X7(xiX8 + X3X6 + X4X5 - X2X7) - 2X3X8X5X7 

X221 = a;g(xiX7 - X3X5) + X7(X6X4 - X2X8) 

^222 = X8(xiX8 - X3X6 - X2X7 - X4X5) + 2X7X8X6X4 



Theorem A. 2. A periodic network of relation with n x n fundamental domain is 
bipartite realizable if and only if it is realizable and at any vertex v, the signature 
satisfy 

vfvl + vlv^ + V^vl + 'U4'i;5 - 2viVsV2V'j — 2viV8V3Ve — 2'i;i'U8'i;4'i^5 

-2v2VrV3VG - 2v2VyV4V5 - 2v3VeV4V5 + AviViVQV-j + 4'i;2V3'i;5'i;8 = 
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Proof. Without loss of generality, we assume at each matchgate, the signature {111} 
is 0, and assume f is a black vertex. The {111} entry of the signature of a black 
vertex is gives us 

aiiai2ai3f 1 + auai2V2 + anflis'^s + anv^ + ai2ai3V5 + ai2VQ + ai^vj + vg = (42) 
The parity constraint implies that the {110} entry is also 0, namely 

011012003^1 + aiiai2f 2 + aiiao3t^3 + «ii^^4 + 012003^^5 + «i2^^6 + aosvr + vg = (43) 

From (142|) ( 143|) . we can solve ail, and solution has following form 

^Ni^N2 
~ D,~ D2 

Then N1D2 — N2D1 = 0. Under the assumption that 093 7^ 013, we have 

{-V2V5 + V6Vi)al2 + {VsVl + VqVs - ViV5 - V2Vj)ai2 + ^8^3 " V^Vj = (44) 

From (12), we have 

2{viVq - V2V5)ao2ai2 + (vivs + V3VQ - - V2V7){ao2 + aia) + 2{v3Vs - v^Vj) = 0(45) 
2x fHil) — (l45l). under the assumption that ao2 7^ ai2, we have 

ViVs + V3VG - - V2V7 



ai2 



2(f - V2V5) 



which implies that equation f l44p has double real roots, and its discriminant is 0. 
That is exactly the statement of the theorem. □ 

Proof of Proposition 3.9 Since holographic reduction is an invertible process, 
two matchgrids M and M are holographically equivalent if and only if there exists a 
basis for each edge, such that one can be transformed to the other using holographic 
reduction. Obviously holographic equivalent matchgrids have the same dimer parti- 
tion function. Assume weights mfe,m^ of M and mfe,m^ of M are as follows: 



( 4' b{' a\' d\' ) 



m;^' = ( 4' bi^ 4' )* 

t 



( ci^' dY d\' ) 



ml 

= ( b^ d^ d^ y 
Then by equations (fT6|) - (!2T]) . and the uniqueness of the basis on each edge, we have 

aoWi3 = boi~^b'{i~^ (46) 

44 = bl2''b\-2''^ (47) 

\ = 44 (48) 
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-1 


44 
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b'ld'l 


i— 

^2 




\)lcl 
a(d^l 


C?2 ^2 
'-'2 ^2 


ij i 

aQ^a 



Plugging in (!46])-(!4HD to ©-([lO!), we have 



t'2 ~ '-2 "'01 "'02 Pl3 '^l 
lij lij i.jfi idfi *iJiO ^ij 

^ij ij i,j,0 i,j,0 i,j,0 y^ij 

0-2 —'^2^11 '"02 '"03 ' ^1 
jij lij i.j.O i.j.O i.j.O ^ij 
«2 = "2 Pll P12 Pl3 ■ Ci 



= 




^jjil ^ijA i^j^^ 
"'01 "02 Pl3 




"-01 2^12 "-03 

«i 


a{ = 
dY = 


i-tJA 
Pll "-02 ''^OS 


i.j.l i.j.l i.j.l 
Pll P12 Pl3 



To prove that probabihty measures of M and M are identical, we only need to prove 
that for any dimer configuration, the products of weights differ by the same constant 
factor. Each dimer configuration corresponds to a binary sequence of length N , 
where = 3?t,^ is the number of connecting edges. Choose an arbitrary edge with 
riQ Til 



basis ( 1. If the edge is occupied, then from M to M, adiacent generator 

Po Pi J 

weight is divided by pi, while adjacent recognizer weight is multiplied by pi. If the 
edge is unoccupied, then from M to M, adjacent generator weight is divided by no, 
while the adjacent recognizer weight is multiplied by hq. Therefore, the total effect 
is for any particular dimer configuration w 



Partition (cu on M) 



n c'r^2^' 



Partition(cu on M) 
which is a constant independent of configuration w. □ 
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